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ABSTRACT 

We analyse the stellar content of almost 300 deg^ of the sky close to the Galac- 
tic plane by directly comparing the predictions of the SKY model (Cohen 1994,1995; 
Wainscoat et al. 1992) with star counts taken from the Two Micron Galactic plane 
Survey (TMGS: Garzon et al. 1993). Through these comparisons we can examine dis- 
crepancies between counts and model and thereby elicit an understanding of Galactic 
structure. Over the vast majority of areas in which we have compared the TMGS 
data with SKY's predictions we find very good accord; so good that we are able to 
remove the disc source counts to highlight structure in the plane. The exponential 
disc is usually dominant, but by relying on SKY's predicted disc counts we have been 
able to probe the molecular ring, spiral arms, and parts of the bulge. The latter is 
clearly triaxial. We recognize a number of off-plane dust clouds not readily included 
in models. However, we find that, whilst the simple exponential extinction function 
works well in the outer Galaxy, closer than about 4 kpc to the Centre the extinction 
drops dramatically. We also examine the shape of the luminosity function of the bulge 
and argue that the cores of all spiral arms we have observed contain a significant pop- 
ulation of supergiants that provides an excess of bright source counts over those of a 
simple model of the arms. Analysis of one relatively isolated cut through an arm near 
longitude 65 degrees categorically precludes any possibility of a sech^z stellar density 
function for the disc. 
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1 INTRODUCTION 

Although star counts have long been used for the deter- 
mination of Galactic structure (see Paul 1993), there are 
still a large number of parameters which are not well deter- 
mined. Reviews of star count surveys can be found in Bah- 
call (1986), Majewski (1993) and Reid (1993). Many visible 
studies have been directed towards the Galactic poles where 
extinction in the Galactic plane can largely be ignored. Star 
counts, velocity distributions and metallicity have been used 
to determine and separate the various Galactic populations, 
such as the old disc, young disc and halo. Detailed reviews 
are given, for example, by Gilmore, Wyse & Kuijken (1989) 
and Freeman (1987), and references therein. Hence, the ver- 
tical structure of the local Galaxy is well measured (although 
there remain contentious issues such as the thick disc). The 
horizontal structures are less well studied. Robin, Creze & 
Mohan (1992) obtained deep visible star counts in the anti- 
centre region and concluded that the disc has a scale length 



of 2.5 kpc and a radial cutoff at about 14 kpc. However, 
apart from individual clear windows, visible star counts in 
the plane are limited to a maximum heliocentric distance of 
about 3 kpc (Schmidt-Kaler 1977), because of the high ex- 
tinction. A method of overcoming this, at least in part, is to 
observe in the infrared, where the extinction is significantly 
lower. Infrared source counts detect sources to a greater dis- 
tance and are far less affected by local dust, so the measured 
distribution of sources is closer to the true distribution. 

To date, infrared surveys of the Galactic plane capable 
of detecting point sources have been limited. The largest 
IR survey to date was that performed by the Infrared As- 
tronomical Satellite, IRAS, which mapped nearly the en- 
tire sky at 12, 25, 60 and 100 fim. Much work on Galactic 
structure has been enabled by IRAS particularly on the 
distribution of very late type stars (see Beichman, 1987, for 
a review). However, at large distances, IRAS could detect 
only extremely luminous sources (e.g. OH/IR stars) and 
was severely confusion-limited over much of the Galactic 
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Table 1. Principal attributes of SKY. For more details see WCVWS, Cohen (1994,1995), and Cohen et al. (1994). 



Geometric 
components 



No. of types Radial scale Source scale Mag. range Wavelength 
of source length height modelled range 



disc; bulge; arms, spurs, Gould's Belt; ring; halo 



87 



3.5 kpc 



90-325 pc (-10,+30) 0.14-35 /im 



Table 2. Summary of regions of the sky analysed. Listed for each strip are the declination and I as the strip crosses the GP. Normally a 
strip is identified by one of these two parameters. Also listed are the starting and ending points for each strip, the number of useful areas 
of comparison with SKY, the area per bin, and the relevant figures in the Appendix for this strip. The area is important as it relates 
observed TMGS counts in an area to counts per deg'^ (and the associated Poisson errors are calculable as shown in Figure 1). 



Dec. 


Plane crossing 


Start 


End 


No. useful 


Typical 


Fig. No. 


(J2000) 
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(;,6) 


(/,6) 


zones 


area (dcg ) 




-30 


-1 


(6.54,-14.5) 


(-5.53,7.5) 


45 


1.9 


A1,B1,C1,D1 


-22 


7 


(15.007,-15.0) 


(-2.223,14.5) 


61 


2.0 


A2,B2,C2,D2 


-11 


21 


(27.495,-15.0) 


(11.619,14.5) 


61 


0.69 


A3,B3,C3,D3 


-5 


27 


(29.294,-4.5) 


(24.124,5.5) 


20 


0.18 


A4,B4,C4,D4 


-1 


31 


(38.994,-15.0) 


(23.264,15.0) 


61 


2.7 


A5,B5,C5,D5 


+3 


37 


(39.74,-5.5) 


(33.968,6.0) 


24 


0.25 


A6,B6,D6 


+19 


53 


(56.215,-5.0) 


(50.582,5.5) 


22 


0.18 


A6,B7,D7 


+27 


65 


(73.688,-5.5) 


(66.817,5.5) 


61 


2.0 


A7,B8,C6,D8 


+32 


70 


(75.837,-15.0) 


(56.995,15.5) 


60 


0.56 


B9,D9 



plane. The Two Micron Sky Survey (Neugebauer & Leighton 
1969: hereafter TMSS) also had a low limiting magnitude 
(rriK ~ 3) and so could not detect sources in the inner 
Galaxy. Other surveys with greater sensitivity have tended 
to be restricted to the Galactic Centre (hereafter GC) re- 
gion. Those that have, at least in part, examined the disc 
(e.g. Eaton et al. 1984; Kawara et al. 1982) have covered only 
very small areas of sky, typically a few hundred arcmin^. 
These surveys ran the risk of sampling an area with un- 
usually high extinction, as well as suffering from low source 
counts, and were unable to show the large-scale distribution 
of sources. 

Garzon et al. (1993) have described the Two Micron 
Galactic Survey (hereafter TMGS). This is a K-h&nd sur- 
vey with an approximate limit of rriK ~10 which represents 
a valuable addition to our knowledge of the near-infrared 
Galaxy, deepening the broad sweep coverage of the plane by 
a factor of almost 1000 compared with the TMSS. It ful- 
fils an important transitional role, providing a homogeneous 
data base over a number of regions in the Galactic plane at 
a time when both the hemispheric DENIS (Epchtein 1998) 
and all-sky 2MASS (Skrutskie 1998) three-band surveys are 
in routine operation. These new surveys will go about 4-5 
magnitudes fainter than TMGS and offer simultaneous ob- 
servations in UK' or JHK' . Nevertheless, a 7^-band survey 
with a magnitude limit of rriK ~10 such as the TMGS can 
usefully probe the structure of our Galaxy. 

Hammersley et al. (1994) and Calbet et al. (1995) ex- 
plored the morphology of star counts from the TMGS by 
heuristic methods, while Hammersley et al. (1995) deduced 
the displacement of the Sun from the plane from the slight 
north-south asymmetry of counts crossing the plane. How- 
ever, these papers are based principally in looking at the 
very obvious features seen in the the TMGS. In order to 
understand the K star counts more fully, a detailed model 
is required which allows the total counts to be broken down 



into counts from the various Galactic geometric components 
in each magnitude range. The model that we have chosen 
to work with is described by Cohen (1994), who has en- 
hanced the "SKY" model for the point source sky originally 
described by Wainscoat et al. (1992: hereafter WCVWS), 
adding molecular clouds for mid-infrared realism, extend- 
ing the Galactic structure to model the far-ultraviolet sky, 
and displacing the Sun 15 pc north of the plane (Cohen 
1995). 

The purpose of this paper is to describe our use of the 
SKY model to provide a formal basis for interpreting the 
TMGS source counts. Our eventual objective is to refine 
the model, guided by any discrepancies that we find with 
the TMGS data base. To do this we contrast the TMGS 
data from almost 600 patches of sky (some 296 deg^) in the 
vicinity of the plane with the predictions of the SKY model, 
to investigate the stellar content and probable geometry of 
the Galactic plane. The advantage of this approach is that 
we can emphasize otherwise slight deviations from expecta- 
tion, and thereby define those regions in which the TMGS 
has uncovered substantive evidence of structure contrary to 
that in a formal model. We present analyses with SKY of 
TMGS star counts along entire survey strips, from 10° to 
30° in length, cutting the plane at a variety of longitudes. 
We conclude that SKY provides a very good overall descrip- 
tion of K counts in the plane, which is chiefly dominated 
by the exponential disc, but that the TMGS offers impor- 
tant insights into the character of structure very close to the 
Galactic plane and the nature of the bulge. We find four ba- 
sic styles of statistically significant deviation of SKY's pre- 
dictions from the observed counts. One appears to arise from 
a significant constituent of the ring, or arm, populations not 
currently in SKY, and strongly confined in latitude close to 
the plane in the inner Galaxy (longitudes 20°-31°). A sec- 
ond type of behaviour suggests that the luminosity function 
(hereafter LF) of the inner bulge is not that implicitly rep- 
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Table 3. Summary of lines-of-sight covered. Listed for each strip is the closest approach to the GC (assuming R0=85OO pc) and the 
features that arc sampled. The question marks with the ring/bar depend on the interpretation of the counts (see Sections 8,9). The Arm 
crossings listed are those closest to the sun (Sct=Scutum, Sgr=Sagittarius, and Per=Pcrseus arm). The * indicates that the line of sight 
is almost tangential to the spiral arm. 
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resented by SKY. The third set of departures from an ideal 
model is caused by spiral arms that do not lie in the plane 
(cf. WCVWS; Cohen 1994). The remaining type of pattern 
is due both to off-plane dark clouds not included in SKY's 
smooth analytic expression of the extinction, and lines of 
sight through the Galaxy which have unusually high or low 
distributed (rather than localized) extinction. 



2 THE TMGS AND THE SKY MODEL 

2.1 The TMGS 

The TMGS (Garzon ct al. 1993) is a collaborative project 
between the Instituto de Astrofisica de Canarias, Tenerife, 
and Imperial College, London, to map large parts of the 
Galactic bulge and plane. Observations were made between 
1988 and 1995 on the 1.5-m Carlos Sanchez Telescope of 
the Observatorio del Teide, Tenerife. Some 350 deg^ of sky 
have been mapped down to about 10 mag in K, resulting 
in a data set of about 700,000 sources. The main areas are 
in strips about 30° in RA (centred on the Galactic plane) 
by about 1.5° in declination. Most TMGS sources lie within 
a few degrees of the Galactic plane, the majority having 
no optical counterparts to at least mv~lG, implying that 
they are sources deep within the Galaxy. The star counts 
presented in this paper were made using all of the data in 
each area. Although the data were taken in strips in RA, the 
TMGS catalogue has been converted to Galactic coordinates 
and the counts made using 1° bins in latitude. 



2.2 SKY: the Model 

The original SKY model described in WCVWS incorporates 
six fully detailed Galactic components: an exponential disc; 
a Bahcall-type bulge; a de Vaucouleurs halo; a 4-armed spi- 
ral with 2 local spurs; Gould's Belt and a circular molecular 
ring. Each geometrical element may contain up to 87 source 
categories, derived from detailed analyses of the content of 
the IRAS sky (Walker et al. 1989). These comprise 33 "nor- 
mal" stellar types; 42 types of AGE star, both oxygen- and 
carbon-rich; 6 types of object that are distinct from others 
only in their MIR high luminosity; and 6 types of exotica 
including H II regions, planetaries, and reflection nebulae. 



Every category has its own set of absolute magnitudes and 
dispersion in ten hardwired passbands, and an individual 
scale height and volume density in the local solar neigh- 
bourhood. Some sources are absent from some components 
(the arms and ring are made rich in high-mass stars; the 
halo is deficient in these). 

We have used version 4 of SKY as detailed by Co- 
hen (1994,1995) and Cohen et al. (1994), the principal at- 
tributes are summarized in Table 1. This is different from 
WCVWS in that it has the solar offset from the plane and 
reduced halo:disc ratio, however, we have not included the 
additional molecular cloud populations essential to obtain 
a high-fidelity representation of the mid-infrared sky. This 
omission is justified because the clouds were inserted by Co- 
hen (1994) for IRAS wavelengths of 12 and 25 fim, which 
are essentially insensitive to extinction. Therefore only the 
extra sources were included and not the extra extinction 
through the clouds, yet the extra populations associated 
with these clouds are surely accompanied by extra extinction 
not accounted for within SKY. Cohen et al. (1994) demon- 
strated in the far-ultraviolet that one could assign a value to 
the additional extinction encountered through local molec- 
ular clouds using SKY. Within the TMGS archive one can 
clearly recognize obscuring clouds beyond which no 2-pm 
sources are detected, although at K only the dense cores 
cause a significant effect on the counts and so the effect of 
the extinction is generally limited to single TMGS zones. 
Consequently, without a detailed prior knowledge of the ad- 
ditional extinction along the specific TMGS lines of sight, 
it is preferable to model without including these molecular 
complexes at wavelengths susceptible to appreciable atten- 
uation through clouds of Av>10 mag. In this way one can 
more readily identify those lines of sight that intercept dark 
clouds at K. Furthermore, the IRAS sources detected in the 
clouds had warm dust which makes the sources very lumi- 
nous at 12 ^m. This warm dust does not have as large an 
effect on the absolute magnitudes at 2.2 pm, so the sources 
will be detected at K but will not be that much brighter 
than the sources without warm dust. Hence they will be 
more difficult to see against the background of disc or bulge 
sources. This is the sole modification made to the published 
description of SKY4. 
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2.3 Why Confront SKY and the TMGS? 

The SKY model was originally built to reproduce the IRAS 
12- and 25-pim source counts. However, it has been suc- 
cessfully compared with observations at many other wave- 
lengths including the far-ultraviolet (Cohen 1994; Cohen et 
al. 1994). The model has already been compared with DE- 
NIS K star counts in small regions (Ruphy et al. 1997) for 
the magnitude range mK=10-13. However, the comparison 
with the TMGS is possibly the most challenging for the 
model because its magnitude range, mx=5-9, probes the 
brightest sources in the LPs of the geometric components 
such as the arms and bulge. At rriK, these LPs are largely 
unknown, particularly for the arms, but the contrast with 
the disc will be maximal so any discrepancy with the model 
will readily be seen. Further, the areas to be studied are long 
cuts through the arms and disc, providing nearby reference 
regions out of the arms. This, too, allows differences to be 
explored in a manner that is impossible given only a sin- 
gle region. The final important advantage is the relatively 
large area of sky covered which, in many regions, affords 
meaningful star counts to be produced as bright as mK=5. 
The fact that SKY produces good matches to IRAS 
(Cohen 1995) and visible stars counts towards the NGP 
(WCVWS) does not guarantee that it will successfully 
match the TMGS star counts. The K LP is dominated by 
photospheres whereas, for IRAS, the principal populations 
detected at significant distances are all dust-enshrouded. 
Clearly, the IRAS sources will be detected by the TMGS 
but are many magnitudes fainter, greatly diminishing their 
importance. Thus, whilst SKY may get the LP of the bulge 
correct for IRAS, this does not necessarily imply correct- 
ness for the TMGS. Garzon et al. (1993) already noted that 
the Galaxy at 2.2 /im looks very different from that in the 
visible because of extinction in the plane. Purthermore, the 
principal sources detected by the TMGS are K4-M0 giants 
or M-supergiants, which are far less important in visible star 
counts when coupled with extinction. 



3 CONFRONTATION OF THE MODEL BY 
THE OBSERVATIONS 

In the Appendix (Figures A1-A7, D1-D9), we present 
TMGS counts along nine strips in the first quadrant. These 
illustrate our interpretations of the counts by direct com- 
parisons with the predictions of SKY. Differential counts 
were scaled per square degree per magnitude, using entirely 
independent TMGS data in 0.2-mag wide bins. 

Table 2 summarizes for each strip: the J2000 declina- 
tion; the longitude of the Galactic plane crossing; the bound- 
ing coordinates; the number of useful (i.e. with tolerable 
implicit Poisson errors) zones we present in the relevant fig- 
ure; the typical area sampled at each position along the strip 
(this varies slightly along a strip); and the figure which illus- 
trates the star counts and SKY predictions. Table 3 amplifies 
this information by indicating the closest approach of each 
line of sight to the Galactic Centre and which (non-disc) 
geometric components contribute to the predicted counts. 

Our analyses depend on three basic analytical tools: dif- 
ferential star counts (Figures A1-A7, D1-D9); profiles of cu- 
mulative counts in a specific AT-magnitude bin as a function 



of latitude along each strip (Figures B1-B9); and similar 
profiles created after subtracting the disc component pre- 
dicted by SKY from the observed counts (Figures C1-C9). 
In both Figures B and C, we represent the observations by 
solid lines and the model counts by dashed lines. 

Figures l(a)-(h) present eight representative fields, 
taken from the montages shown in the Appendix, for which 
we explicitly plot the Poisson error bars appropriate to the 
original counts (which were scaled by 5, to yield counts per 
mag from the 0.2-mag binned data, and again, according to 
the actual area sampled, to yield counts per deg^). Figure 
1 distinguishes the counts arising from separate geometric 
components as follows: solid line - total count; dotted or less 
heavy line - disc; long dashes - spiral arms, local spurs and 
Gould's Belt; short dashes - bulge; short dashes and dots 
- molecular ring; long dashes and dots - halo. The figures 
were chosen to be representative of the various regions cov- 
ered as well as showing specific features in their form which 
are discussed later in the paper. 

As with all surveys, Malmquist bias can arise in TMGS 
and we see its effects in Figure 1(h) where the counts in the 
final few faint bins rise abruptly before the limit of complete- 
ness is attained. Note also that this is a region in which the 
actual source density is relatively low. In general, the direct 
comparison of SKY and TMGS counts indicates that TMGS 
is complete to a magnitude between 9.2 and 9.8, depending 
on the local surface density of sources. 



4 EXTINCTION 

4.1 Off-Plane Dark Clouds 

As with all previous versions of the model, SKY4 represents 
interstellar extinction simply by the product of two expo- 
nentials; one radial, with scale length of 3.5 kpc (matching 
that of the stars); one vertical, with scale height of 100 pc 
(matching that of the youngest stars still associated with 
dark clouds) . This version of SKY does not include off-plane 
dark clouds so, if a TMGS strip runs through the core of 
such a cloud, the counts will be significantly decreased when 
compared with the model (however, the greatly reduced ex- 
tinction at K compared with the visible implies it requires 
many magnitudes of visible extinction to cause a significant 
effect on the TMGS counts). Figures A5,B5,D5 illustrates 
this phenomenon for the area centred on I = 29.3°, 6 = 4° 
from the I = 31° strip. One sees that the TMGS counts fall 
well below the predictions of SKY yet, at & = 2.5° or fe = 5°, 
the model and the TMGS are in good agreement. Figure 
1(a) enlarges the plot for 6 = 4°, the deepest discrepancy, 
and adds the Poisson error bars. This location corresponds 
to a known off-plane cloud which clearly shows up in emis- 
sion in longer-wavelength COBE maps, whereas the area is 
almost devoid of visible stars. In fact, the core of this cloud 
does not fill the TMGS bin at this location and sampling 
at higher resolution would decrease the TMGS counts even 
more at the centre of the cloud. The only potential H I coun- 
terpart (Weaver 1998, priv. com.) is a fairly uniformly dark 
patch (i.e. an absence of H I) well-centred in an H I emission 
ring, at an LSR velocity of —15 km s^^, and about 1-1.5° 
in diameter. Figure 1(b) shows another example, this time 
of a more localized off-plane cloud, whose influence is recog- 
nizable by an abrupt change in the observed source counts 
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Figure 1. g-h. Differential star counts versus magnitude sliowing the detailed comparisons between the TMGS and SKY for sample 
regions. The black diamonds are the TMGS points with Poisson error bars marked. The SKY geometric components arc marked as 
follows: solid line — total count; dotted or less heavy line - disc; long dashes - spiral arms, local spurs and Gould's Belt; short dashes — 
bulge; short dashes and dots - molecular ring; long dashes and dots — halo 



fainter than a specific magnitude. Here one sees an offset in 
source counts by 0.2 dex for rriK >8.5. Analysis of this fine 
of sight by SKY indicates that the dominant contributors at 
this magnitude are disc KO-3 giants, at a probable distance 
of about 700 pc. Consequently, we suggest that one encoun- 
ters a cloud with incremental Ay ~ 5 mag in this direction 
and at this distance. Again, H I maps (Weaver 1998, priv. 
comm.) show a sizeable void, bigger than that seen above 
(about 2° in diameter), well-centred on this TMGS field, in 
an emission ring at an LSR velocity of -1-38 km/s. In both 
cases, the presence of H I emission with central voids sug- 
gests that the hydrogen in these clouds is largely molecular 
rather than atomic. 



4.2 On-Plane Dark Clouds and General 
Extinction in the Plane 

The TMGS reveals a number of "on-plane" (|fe| < 2°) dark 
clouds, particularly towards the inner Galaxy. Hammersley 
et al. (1996) discussed the dark clouds which cause the major 
loss of stars at I = —1°, b — 0°, and showed that most of 
the extinction comes from individual dark clouds close to 
the GC. However, there are individual clouds on, or very 
near, the plane in many TMGS areas which locally cause a 
noticeable loss of counts (/ = 7°, 6 = 1° or / = 31°, 6=1°: 
Figures A2, B2, D2, A5, B5 and B5) 

SKY includes only the general extinction along the line 
of sight. This is modeled as a double exponential which con- 
tinues all the way into the GC but is normalized in the 
solar neighbourhood. Figures Bl- B9 show the cumulative 
TMGS and SKY star counts (the latter separated into to- 
tal and disc contributions) for various strips, as a function 
of Galactic latitude, and the ratios of these totals at differ- 
ent limiting magnitudes. The effect of extinction on SKY's 
counts is clearly visible in these cuts as a dip in counts at 
6 = 0°. For lines of sight with I > 30° , the dip on the plane 



is small, a few percent of the total counts, and there is very 
good agreement with the TMGS star counts. In the inner 
Galaxy, however, SKY predicts a major loss of counts on 
the plane, principally due to extinction in front of the bulge 
sources (the difference between SKY's total and disc curves 
for rriK > 7 and the I — 7° and —1° strips is almost to- 
tally due to the bulge). The model does reasonably well near 
/ = — 1°, 6 = 0° although, as already stated, most of the ex- 
tinction appears to come from identifiable dark clouds close 
to the GC, not represented in SKY. However, SKY substan- 
tially over- predicts the extinction on the plane at / — 7°, 
6 = 0°. The ratios of model-to-observed star counts for this 
region show a narrow peak which the original plot indicates 
is due to extra extinction in SKY. 

Hammersley et al. (1998) show that the extinction at 
the radius of the molecular ring {OARq < R < 0.5Rq) 
is about a factor of three above that predicted by a sim- 
ple exponential model yet SKY's analytical function for the 
extinction in the outer Galaxy (e.g. R> O.5R0) gives a rea- 
sonable representation for the distributed extinction in the 
plane. Therefore, the only explanation for why SKY so sig- 
nificantly overestimates the extinction in the line of sight 
towards / = 7°, 6 = 0°, is that the extinction for R > 0.4i?Q 
drops to almost zero until a few hundred pc from the GC, 
where the features associated with the Centre begin (e.g. the 
molecular disc). The idea of there being a hole in the extinc- 
tion in the inner Galaxy is not new: the CO maps clearly 
show a significant fall-off in density inside the molecular ring 
(Clemens, Sanders & Scoville 1988; Combes 1991). Freuden- 
reich (1998) also needed a hole in the extinction in the inner 
Galaxy to match the DIRBE surface brightness maps. 



5 THE EXPONENTIAL DISC FOR B > 5° 

SKY incorporates a single exponential disc with 87 cate- 
gories of point source, each with its own scale height, colours. 
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and local space density. The scale length of all sources is set 
to 3.5 kpc. SKY does not include a thick disc for several 
reasons. Firstly, it provides excellent fits to counts at a wide 
variety of wavelengths without using a thick disc. Secondly, 
the proportion of disc counts contributed by the thick disc, 
even at optical wavelengths where metallicities and kine- 
matics may help separate these two populations, is minimal 
(< 2%; e.g. Gilmore, Wyse & Kuijken 1989). Thirdly, the 
stars of the thick disc contribute negligibly in the infrared, 
far below the levels of Poisson errors of the TMGS (and 
many other) data sets. 

SKY separates young and old components by assigning 
the young population (OBA dwarfs and supergiants) a sig- 
nificantly smaller scale height, 90 pc, than the 325 pc of the 
old populations. Hence, the vast majority of the sources that 
SKY predicts that the TMGS should see in the exponential 
disc are the old K- and M-giants for regions well away from 
the bulge and more than 5° from the plane. The relevant 
figures are those for the strips ai I = 21°, 31° and 65° (Fig- 
ures A3, A5, A7, D3, D5 and D8). Figures B3, B5 and B8 
show the associated cumulative counts for the TMGS and 
the model for these regions, to tuk = 9. 

Both the form of the curves and the differential star 
counts show that SKY gives a good fit to the TMGS at all 
magnitudes in these longitudes, for |&| > 5°. 

In most areas of the sky, the disc dominates the star 
counts, making the remaining geometric components more 
difficult to identify and analyse. However, the accuracy of 
the fit of SKY'S disc to the TMGS counts makes it en- 
tirely reasonable to subtract the model's disc counts from 
the TMGS counts. This highlights the distribution of sources 
in the arms, ring, and bulge, if accomplished with care. 

Alternatively, some argue that an exponential disc is no 
longer justifiable at the lowest latitudes in the Galaxy, and 
that the sech^z form might be physically more defensible. 
There are a number of papers comparing the seda^z and ex- 
ponential in the discs of external galaxies (e.g. van de Kruit 
1988; De Grijs, Peletier & van de Kruit 1997). However, the 
results, even in the infrared, are complicated by extinction 
when close to the plane of the galaxy. The sech^z function 
asymptotically approaches an exponential as stellar densi- 
ties fall far from the plane but, for latitudes within about 
10° of the plane, it fails utterly to reproduce the observed 
counts. We illustrate this failure in Figure 2. The 1—65° 
strip is one in which we see only disc and arms, the lat- 
ter falling almost two orders of magnitude below the counts 
from the former. Thus, it is an ideal location to experiment 
with a sech^z distribution function. Figure 2 shows the cut 
through the plane at m_K-=9 for two model discs for /=65°, 
\b\ < 15°, one is exponential in height (this is a simplified 
model based on SKY but it is not the SKY model) and the 
second is identical except that it has a sech'^z height distri- 
bution. We have normalized the sech^z curve to the wings of 
the cumulative latitude profile in the regions with |&j > 10° 
as both functions work well, away from the plane. Yet, for 
all inner latitudes, the sech^z function delivers no more than 
half the requisite counts. This is equivalent to saying that, 
everywhere in the disc, even in those locations where spiral 
arms and other components do not make appreciable con- 
tributions to counts, more than 50% of the observed stars 
cannot come from the disc. Were we to have carried out this 
test at other TMGS positions in the plane, the excess of 
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Figure 2. Two simplified (not SKY) model euts at (=65° for 
mx=9. The solid line shows an exponential disc (as in SKY) the 
dashed is based on a sech^z. In both cases the SKY disc LF was 
used. The two curves have been normalized at 6=10°. 



counts above the sech^z law would be even larger than that 
shown in Figure 2. 

We further note that the SKY model produces good 
fits to the DENIS K star counts in the plane to mx=13.5 
(Ruphy et al. 1997), significantly deeper than the TMGS. In 
regions away from the Galactic bulge, the contribution from 
the disc will be at least an order of magnitude greater than 
from any other component, e.g. the arms. If the disc were 
not exponential, or at least very close to exponential, all the 
way to the plane, then SKY would significantly overestimate 
the counts on the plane. 

The sech^z function is strictly valid only for an isolated, 
self-gravitating, isothermal system. Self-consistent dynami- 
cal models of the disc (e.g. Haywood et al. 1997) show the 
dependence of the resulting z-distribution on star formation 
rate (SFR). For all but declining SFRs, the vertical distri- 
butions can mimic exponentials. It is our suspicion that the 
sech?z function from a locally isothermal disc achieved its 
popularity in the context of overall light distributions across 
the planes of edge-on spirals at optical wavelengths (e.g. van 
der Kruit & Searle 1982). In such a situation, heavy and 
patchy extinction militates against any meaningful compar- 
ison between a sech^z law and observations and, far from 
the planes of such galaxies, one cannot distinguish between 
sech^z and exponential functions. Perhaps the only redeem- 
ing feature of a sech^z law is that it offers a null second spa- 
tial derivative at zero latitude, yet this is precisely the lati- 
tude range in which optical observations are invalid for the 
purposes of making such comparisons. Further, even in the 
field of edge-on external galaxies, there is evidence that the 
infrared light distributions follow exponential rather than 
sech^z functions (e.g. Wainscoat et al 1989). 
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6 THE PLANE OF THE DISC 

For |fe| < 2°, at all longitudes, the TMGS has more counts 
on the plane than are predicted by SKY. The profiles of 
counts (Figure B) show three distinct distributions at these 
latitudes: 

• A distribution extending from about |fej < 4° that is 
seen at all longitudes between / = 70° and —1°. 

• A spiked distribution about 1° wide centred on the 
plane at I = 21° and 27°. 

• The bulge at I — —1° and / — 7° (which we treat in 
section 9). 



6.1 The Spiral Arms 

SKY predicts that the 2-^m star counts should form a con- 
tinuous distribution from b — 15° to 6 = 1° with the only 
serious deviations due to extinction in the plane. However, 
the scans I = 21°, 31° and 65° all show a distinct transition 
at about |fe| = 5°. (These can be seen both in the cuts [Fig- 
ures B3, B5 and B8] and in the counts [Figures A3, A5, A7, 
D3, D5, D8 ].) The star counts for |fe| > 5° approximate an 
exponential very well but, at fe = 5°, it is clear that a sec- 
ond contribution is required with a far smaller scale height 
than typifies the stars of the disc (i.e. an element rising more 
rapidly than the disc as latitude decreases) . Figure 1 (e) rep- 
resents this phenomenon from the 1—65° strip, some 3.5° off 
the plane. 

A feature so close to the plane has to be due to a young 
component. There are two possibilities to explain this dis- 
tribution; a young disc or spiral arms. SKY already has a 
young component within its representation of the disc and 
these young stars (O, B and A) have a much smaller scale 
height than the older stars. Thus, if the extra sources in the 
plane were young disc stars, this would imply an enhance- 
ment in these populations. Indeed, examination of Figures B 
and C shows that, particularly at the brighter magnitudes, 
the numbers of young sources would need to be increased by 
orders of magnitude to reproduce the counts near the plane. 
This would create a major distortion of the LF, clearly rec- 
ognizable in the local vicinity. Figure C shows the TMGS 
counts after subtracting the SKY disc. Note that, for both 
I — 31° and I = 65°, the width of the non-disc source distri- 
bution is independent of magnitude, suggesting that these 
sources are grouped at a specific distance rather than spread 
along the line-of-sight. For the above reasons, we prefer to 
attribute the excess of stars on the plane to the arms, rather 
than to a boost in the young disc. 

The prediction from SKY is that the arms are always 
substantially below the disc counts, although there are cer- 
tainly locations in which the arms' contribution is second 
only to that of the disc. However, particularly at the brighter 
rriK magnitudes {rriK < 7.5) the TMGS finds that the counts 
from the arms are at least as important as those from the 
disc within 2° of the plane. The two major TMGS regions 
(i.e. with sufficient area to give reliable statistics at the 
brightest magnitudes) that run only through spiral arms and 
through none of the other inner Galaxy components are at 
I = 31° and I = 65°. At I = 31°, 6 = 0° at niK^T, the 
TMGS has nearly three times SKY's predicted counts. At 
/ = 65° at mK~7, the excess on the plane is only 50% but 



at mK~5 the TMGS is again a factor of 3 too high. This 
can be clearly seen in Figure B. 

By mK=Q in both the I = 31° and 65° strips the arm-to- 
disc ratio has dropped significantly. In the range 8 < 7^ < 9 
at & = 0° there are some 331 predicted disc sources compared 
with only 40 arm sources. At I = 31° the loss of contrast is 
not as great (540 disc sources to 350 arm sources) probably 
because this line of sight is almost tangential to the Scutum 
arm, looking into the inner Galaxy, whereas the I — 65° line 
of sight roughly maintains a Galactocentric distance of Rq 
for over 8 kpc. 

Another interesting conclusion can be reached by de- 
tailed comparison of the star counts at I — 31° and 65°. In 
Figures C5 and C6, the dominant contribution to the disc- 
subtracted counts is from the arms. At niK ~ 5, we see 
roughly the same absolute number (within a factor of 1.5) 
of luminous non-disc sources in both directions yet, by the 
fainter magnitudes (e.g. rriK = 9), the / — 31° line of sight 
shows about 6 times as many spiral arm sources as towards 
/ — 65° . Thus, the LF of the spiral arms must be a function 
of Galactocentric distance, rather than constant as in SKY4. 
If this is truly the case, one might suppose that metallicity 
was implicated. There is a precedent for this argument in 
the observed confinement of the latest WC-type Wolf-Rayet 
stars to the immediate environs of the GC, a phenomenon 
attributed by Smith & Maeder (1991) to the high metallic- 
ity near the Centre which biases evolution of massive stars 
to the formation of WC9 stars. 

At I = 65° (Figure C6), only the Perseus arm is crossed 
and at a Galactocentric distance of about 8.5 kpc. Here the 
excess is around 80 sources deg~^ to mx=8. It is noticeable 
that the peak is not at fe = 0°; rather the centre of the 
distribution is shifted some 1.5° to positive latitude. Also, 
the observed counts drawn from the / = 37° strip are parallel 
to, but above, SKY's predictions at all magnitudes. Both 
features are symptomatic of directions in which the spiral 
arms do not lie in the formal Galactic plane but arch above 
it. This same effect was detected in the IRAS counts by 
WCVWS, who were able to match it to specific regions in 
the H I sky (e.g. Weaver 1974) where neutral gas lies almost 
1 kpc off the plane. They cited, in particular, the I — 50° to 
90° region, which includes the off-plane peak of TMGS K 
counts along the I = 65° strip. 

Figure 3 shows the peak source count detected in each 
strip after the SKY disc has been subtracted. Vertical 
dashed lines are the points tangential to the Scutum and 
Sagittarius spiral arms. The aim is to show primarily the 
arm sources on the plane, therefore, at / = 7° and I = —1° 
only sources to tuk ~ 6 have been shown to avoid bulge 
sources. Also because of the spike (Section 6.2) the values at 
I — 21° have been estimated as if the spike were not there 
and the region at I = 27° has been ignored. The I = 31° 
and I = 37° lines of sight are very close together, yet the 
/ = 31° region has almost twice the number of arm sources 
as at / = 37°. The only significant difference is that the 
/ = 31° direction crosses the Sagittarius and Scutum arms 
whilst / = 37° intercepts only the Sagittarius arm. The ex- 
tra sources to mx=7 and 8 suggest that just under half the 
excess sources can be attributed to the Scutum arm (i.e. 
about 350 sources deg~'^ to mx=8 on the plane). The effect 
of the arms on the counts can be clearly seen in all regions, 
in particular the influence of the number of arms which cut 



© 0000 RAS, MNRAS 000, 000-000 



TMGS vs SKY 9 




80 



40 



20 




Y iK = + 6 \ 
^ AK=+5 I 



60 



40 20 

Galactic LonsiLude 



Figure 3. The peak cumulative count at eacli TMGS plane 
crossing after subtracting the SKY disc. The aim is to show those 
counts from the arms therefore the regions at giving a major 
contribution at fainter magnitudes. Also at 1=21° the counts are 
the estimated counts if the "spike" (see below and Figure 3 of 
Hammersley et al. 1994) arc ignored. 



the line of sight. Directions crossing two arms have nearly 
twice the counts of those crossing only one. Furthermore, 
the inner spiral arms are far stronger than the arms at the 
distance of the solar circle (cf. the difference between Z=65° 
and the other regions) . Another curious feature of Figure 3 
is that at mK=5 and 6 the number of spiral arm sources 
does not alter significantly with longitude between Z=— 1° 
and /=31°. Clearly there is a strong interplay between the 
distance to the arm and the LF but one would expect the 
l=—l° and 1=1° cuts to be significantly less than the Z=31° 
which is almost tangential to the Scutum arm. 

We attribute these excess sources observed near the 
plane by the TMGS to the presence of an additional pop- 
ulation in the arms, not present in SKY4. From the esti- 
mated distances to these spiral arms, this population would 
have to be of supergiant luminosity. Consequently, the spi- 
ral arms clearly provide an important contribution to the K 
star counts, particularly towards the inner Galaxy. 



6.2 The "Spikes" 

The most dramatic deviations from the SKY model any- 
where in the area covered here are the "spikes" on the plane 
at 1=21° and ^=27°. Figures A3 and A4 portray the diflter- 
ential star counts at ^=21°, 6=0°, and ^=27°, &=0°, where 
the excess is substantial at all magnitudes. If the SKY disc is 
removed (Figures C3 and C4), then the TMGS has 6 times 
the counts of the remaining SKY components at /=21°. In 
the cuts (Figures B3 and B4), it can be seen that the distri- 
bution of the excess counts is very narrow, about 1° wide, 
hence the term "spike" . It is significant that the spike only 
appears at 1=21° and 1=27° , there is no evidence for it at 
l=?,l° (or larger longitudes). Neither is there evidence for 



the spike at 1=1° or —1° although the arms are clearly visi- 
ble at the brighter magnitudes in both strips so the absence 
of the spike cannot be just a line of sight effect. This almost 
rules out the feature being due to a continuous feature cir- 
cling the inner Galaxy. The implications of the spikes were 
discussed by Hammersley et al. (1994), who argued that 
these must be features associated with the inner 3-4 kpc of 
the Galactic disc. The scale height must be around 40 pc, 
suggestive of a very young population, and the Mk must at- 
tain at least —11 to —13. Mikami et al. (1982) suggested that 
there is a clustering of supergiants at 1=27° , b=0° . Garzon 
et al. (1997) obtained visible spectra of some of the brightest 
TMGS sources with optical counterparts, proving that there 
are a very high number of supergiants in this region. There 
are two possible explanations for what causes these spikes: 
bar or ring. 

In order to match the IRAS 12- and 25-^im counts 
in this region WCVWS introduced another component, the 
molecular ring (Dame et al. 1987). However, it is clear that, 
even with this representation of the ring, SKY4 is well below 
the TMGS counts. The most significant deviations are con- 
fined to |6|<1.5 ° on the /=21° strip, but persist to higher 
latitudes along the 1=21° strip (|6|<3.5°). Note that the 
slope of the observed counts in these regions more closely 
parallels that of SKY's ring component than the steeply 
sloped disc counts. Were we to posit that an additional stel- 
lar population attends the ring, close to the plane, then this 
population must contribute some 5 times SKY's predicted 
ring counts near 1=21° but almost an order of magnitude 
more than SKY's ring near 1=21°. The excess of sources 
above SKY's expectation is very similar on the 1=21° and 
21° strips, but shows a conspicuous additional bias toward 
very bright sources near 1=27° . 

We note that Ruphy et al. (1997) describe an analysis 
of DENIS near-infrared counts using a new version of SKY, 
with a non-uniform elliptical ring that replaces the old cir- 
cular molecular ring in SKY4. In a later paper we plan to 
exercise SKY5 (after implementing all the changes discussed 
in the present paper) to investigate what TMGS counts tell 
us about the geometry of this ring. In brief, DENIS counts 
probe the lower-density portions of this ring and help to 
constrain the overall geometry of the ellipse. But TMGS of- 
fers direct lines of sight across the first quadrant dense ansa 
included in SKY5. Given sufficient TMGS strips, and new 
deeper JHK images, the density enhancements around this 
new ring should be definable. The inclusion of these star-rich 
ansae provides SKY with a measure of asymmetry (triaxi- 
ality) not present in earlier versions, and can be used to 
provide a more formal examination of the bar. However, we 
defer this discussion to a future paper using new data and 
further developments of SKY. 

Hammersley et al. (1994) suggested that such spikes 
could be associated with a possible bar distribution on the 
grounds that, where the bar interacts with the disc, there is a 
high probability of strong star formation. Note this "bar" is 
not the triaxial bulge, often called a bar, but rather is a stick- 
like bar with a half-length of nearly 4 kpc and an angle of 75° 
to the GC-sun line. Neither can the star formation be related 
to the triaxial bulge somehow triggering star formation in 
the inner disc (as suggested by Freudenreich 1998) as the 
geometrical constraints alone make this nearly impossible. 
Hammersley it et al. (1994) also examined the COBE 2.2- 
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jim surface brightness maps and argued that such a bar far 
more simply and naturally predicts the observed features 
(both images and star counts), whilst the form of a ring 
would have to be contrived using ad hoc solutions to make 
it fit. 

It is beyond the scope of the current paper to engage 
in a discussion of ring versus bar. To make any significant 
advance, new data are required. However, it should be noted 
that many external bars have strong ansae of star formation 
at either ends and so the discussions in the previous two 
paragraphs are not mutually exclusive. 



7 THE BULGE 

Two TMGS strips run through the bulge: those at Z=-l° 
and 1=7° . However, when examining these two regions, one 
must remember that the strips cut the plane at an apprecia- 
ble angle so there is a significant variation in I as well as b, 
which is important this close to the GC These are the most 
complex regions to analyse because disc, arms, and bulge 
contribute to the counts. The first step is to subtract the 
SKY disc to highlight the bulge population. This is entirely 
reasonable because SKY gives a good fit to the exponential 
disc, and because of the high bulge-to-disc ratio fainter than 
mK=7 (e.g. near the plane, there are twice as many bulge 
stars as disc stars). 

Figures Bl and B2 show the cuts for the two strips; Fig- 
ures CI and C2 show the residuals after subtracting SKY's 
disc from the TMGS counts. At ?nx=5 and 6 the residuals 
in the TMGS are clearly due to spiral arm sources because 
they are as narrow in spatial extent, and as large in size, as 
the residuals at 1=31° . However, at mK=6, SKY's predicted 
residuals clearly show a wider distribution than seen in the 
TMGS, due to the predicted bulge sources. Even at 1=7° 
the SKY residual is far wider than seen in the TMGS, yet 
by mK=9 the residuals are far closer together. The explana- 
tion for these discrepancies is that SKY's bulge LF has too 
many of the most luminous sources. A quantitative analysis 
of the K LF, based on the TMGS, is presented in Lopez- 
Corredoira et al. (1997,1999) and shows that SKY's bulge 
LF starts 1 to 2 magnitudes too bright at K. Consequently, 
it is planned to implement a new LF for the bulge and to 
re-examine this issue in the future. 

The SKY bulge is axially symmetric, with a projected 
axial ratio of 0.6. Therefore, if one were to draw ellipses 
of constant predicted bulge star counts, then these would 
have a ratio of major-to-minor axis of 0.6; i.e. all positions 
on the locus VOlG^P+lP would have the same star counts. 
The bulge is the only component which has this type of pro- 
jected spatial distribution and so plotting this function is an 
unambiguous test of whether the TMGS is actually detect- 
ing bulge sources. It should be noted that as the function 
VOlPP+P^ is always positive this means that the line for 
a TMGS strip will double back on itself as it crosses the 
plane. Figure 4 shows SKY's bulge counts plotted against 
y/rat'^P + b^ (right column) together with the TMGS strips 
(left column) through I = 7° (dashed line) and / = — 1° 
(solid line) for various ratios, rat, between 1.5 and 0.2. In 
both cases this is after subtracting the SKY disc and for 
mx=9 so that the bulge will dominate all the remaining 
components. More than a few degrees from the GC, the 



best fit for the bulge (i.e. the lines at the same position 
most nearly overlie each other is for a projected axial ratio 
of 0.6). There is some evidence from the central few degrees 
(Figure 4, l=—l°) that the fit is better for an axial ratio 
of 1.0, but this is consistent with a triaxial bulge (as will 
be discussed below), which would depress the portion from 
negative longitudes. 

Figure 5 shows the same type of plot as Figure 4, this 
time varying the A'-magnitude at a fixed ratio of 0.6. The 
left column shows the TMGS whilst the right column shows 
the ratio of (TMGS-SKY disc) to (SKY total-SKY disc). At 
mK=5, where there is no bulge component, there are two 
separate peaks caused by the spiral arm component in the 
plane. Spiral arms clearly will not follow the bulge locus. 
At m,K=6, SKY (as discussed above), predicts that there 
should be a bulge-like component but the TMGS plots show 
that the only significant population is that of the arms. For 
this reason, apart from the plane crossings, the TMGS is 
below the counts predicted by the model. mK=7 is a transi- 
tional magnitude where the bulge begins to become impor- 
tant, although in the plane there is still a significant arm 
component. On average, the TMGS residuals are below the 
predicted SKY residuals; i.e. the TMGS bulge counts are be- 
low those of SKY. On average at mK=S and 9 away from the 
plane, the ratio of SKY to TMGS is about 1 suggesting that 
by these magnitudes the SKY LF is close to that observed 
by the TMGS. Note that the peak in the ratio where the 
TMGS strips cut the plane (where the curves for each strip 
double back on themselves) at mx=8 and 9 is caused by 
the model overestimating the extinction in the inner Galaxy 
(see Section 4.2). 

SKY predicts that the bulge distribution in the two cuts 
should be asymmetric due to the change in I as well as b. 
However, a comparison of the TMGS residuals, taking ac- 
count of the discrepancy in the LF, shows that for 6 < 0° 
(I > 0° by comparison with the point at 6 = 0°) the TMGS 
is invariably above SKY, whereas for & > 0° (/ < 0°) the 
TMGS is weU below SKY (Figures Bl, B2, CI and C2). It 
should be noted that the effect is very obvious particularly 
at rnK=7 and 8 and so we are not discussing some subtle 
feature in the counts. Initially it was supposed that this was 
due to the extinction by major dark cloud systems above 
the Galactic plane in this direction. However, recent studies 
have shown that in general there is insufficient extinction in 
these clouds to account for the observed effect we see (Ma- 
honey 1998). Typically, the additional Ay above the plane, 
as opposed to below it, is 1 mag (i.e. an extra Ak of about 
0.1 mag). The asymmetry seen in the TMGS would require 
a difference of extinction in K of at least 0.5 mag, which is 
impossible to attribute to extinction alone. Figure 5 shows 
that at l=—l° (the solid line) the difi'erence between the 
TMGS/SKY ratio for & < 0° (Z > 0°, the upper part of the 
solid line) and fe > 0° (Z < 0°) increases with distance from 
the Galactic plane (between 2 and 7 on the x-axis). However, 
the total extinction, as well as the difi'erence in extinction 
above and below the plane, must decrease away from fe = 0°. 
Hence, if the effect were due to extinction both the upper 
and lower parts of the solid line would come together with 
increasing distance from 6 = 0°, and not diverge as is seen. 
Furthermore, if the differential extinction above and below 
the Galactic plane were significant, it would clearly show in 
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Figure 5. The left column shows the TMGS counts after subtracting the SKY disc plotted against {O.d^fi + b^)''-^ for various ratios K 
magnitudes. The right column shows the ratio of (TMGS - SKY disc) to (SKY total - SKY disc ) also plotted against (O.e^P + b2)0.5 
The strips are for l=—l° (solid line) and 1=7° (dashed line). 
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the COBE 2.2-fira surface brightness maps, which is not the 
case (see Freudenreich 1998). 

A closer analysis of the difference between SKY and the 
TMGS shows that, at rriK—S and 9 (where the LFs are in 
closer agreement), the TMGS is above SKY for I > 0° and 
below SKY for I < 0° . Furthermore, at mK=7, where the 
asymmetry in the TMGS is the greatest, the distribution is 
more bulge-like for 6 < 0° (i.e. greater I) than 6 > 0°, where 
the distribution is almost narrow enough to be due solely 
to spiral arms, particularly in the I — —1° strip. A num- 
ber of papers have argued that the bulge is triaxial with 
the near end in the first quadrant at an angle of 12° to 30° 
to the GC-Sun line of sight (Binney et al. 1991: Dwek et 
al. 1995: Freudenreich 1997: Lopez-Corredoira et al. 1997). 
This would on average make the bulge sources at / > 0° 
closer than at I < 0° . If the bulge were triaxial, we would 
expect asymmetries in the TMGS star counts. The sources 
from the closer part of the bulge (/ > 0°) would appear at 
a significantly brighter magnitude than those with I < 0° . 
A rough estimate suggests a magnitude difference as high 
as 0.5 mag. The steep rise in the bulge LF for the bright- 
est sources means the counts cut on very sharply around 
m.K=G.5, but then rise less steeply from rriK >7.5. Thus we 
would expect to see a very asymmetric distribution in the 
mK=7 cut but a far more symmetric one at mx=9, precisely 
as seen in the TMGS counts. Similarly, the ratios shown in 
Figure 5 have the features that one would expect when ra- 
tioing star counts from a triaxial bulge with those from an 
axisymmetric model; i.e. the more positive the longitude, 
the more the counts rise above the predicted values and the 
more negative the longitude, the more the counts fall below 
the predictions. K star counts give access to the LF so we 
suggest that an analysis of the bulge based on star counts is 
likely to be far more sensitive than one based solely on sur- 
face brightness maps because it avoids the ambiguities in the 
size and orientation of the triaxial bulge inherent in an anal- 
ysis of its surface brightness (Zhao 1997). Recently Lopez- 
Corredoira et al. (1997,1999) have analysed the TMGS star 
counts by inverting the fundamental star count equation - 
very different from the qualitative analysis presented here - 
and have found that the bulge is triaxial with an angle to 
the GC-Sun line of 12°. 

For the innermost strip {I — —1°) almost sampling the 
GC, the latitude range in which there is statistically signif- 
icant deviation from SKY that cannot arise from off-plane 
dark clouds, is confined to 1 6|<2.0° (Figures Al, Bl and Dl). 
From these plots, it is apparent that, if the bulge compo- 
nent were simply enhanced by factor of 2, the total predicted 
curve would much more satisfactorily match the counts. This 
bears directly on the bulge LF implicit within SKY. One 
can see from Figure 1(g) how it is possible, quite directly, to 
probe this LF because the observed counts are well-matched 
by SKY except in the restricted range 8.4 > tuk > 6.2. 
TMGS counts exceed SKY's here by about 50% at most, 
but this is naturally attributable to the shape of the bulge 
LF which declines abruptly at bright K magnitudes. We 
note, however, that there is some interplay between changes 
in the density function and changes in the bulge LF. 



8 CONCLUSIONS 

We conclude that, overall, SKY provides a very reliable sim- 
ulator of the TMGS Galaxy. Only rarely are the predictions 
wrong by as much as 50% and, in all these cases, there is 
a fairly obvious change that can be implemented in SKY 
to match the observations by more realistically represent- 
ing the astrophysical situation. Consequently, the wide-area 
coverage offered by TMGS provides substantive constraints 
on the model. This potential interaction between SKY and 
the real Galaxy is extremely valuable in the development of 
a better model and in better understanding the structures 
in the inner Galaxy and their relative importance to K star 
counts. 

Specific changes that we plan to investigate in more de- 
tail and to implement in a new version of SKY are: the inclu- 
sion of an additional population of supergiants in the cores 
of the spiral arms; a modified LF for the bulge; a triaxial 
bulge; more non-coplanar arms and plane than represented 
by WCVWS; and the non-uniformly dense elliptical ring. 

Perhaps the single most important result is that the 
star counts at 2.2 ^m can clearly be modelled using simple 
analytical functions to describe the 4 Galactic components 
of disc, arms, bulge and ring. 



9 



ACKNOWLEDGMENTS 



The Carlos Sanchez Telescope is operated on the island 
of Tenerife by the Instituto de Astrofisica de Canarias in 
the Spanish Observatorio del Teide of the Instituto de As- 
trofisica do Canarias. 



REFERENCES 

Bahcall J. N., 1986, ARA&A, 24, 577 

Beichman C. A., 1987, ARA&A, 25, 521 

Binney J., Gerhard O., Stark A., Bally J., Uchida K., 1991, MN- 
RAS, 252, 210 

Calbet X, Mahoney T., Garzon F., Hammersley P. L., 1995, MN- 
RAS, 276, 301 

Clemens D. P., Sanders D. B., Scoville N. Z., 1988, ApJ, 327, 139 

Cohen M., Sasseen T. P., Bowyer S., 1994 ApJ 427 848 

Cohen M., 1994, AJ, 107, 582 

Cohen M., 1995, ApJ, 444, 874 

Combes F., 1991, ARA&A, 191, 195 

Dame, T. M. et al. 1987, ApJ, 327, 706 

De Grijs R., Peletier R. F., van de Kruit P. C, 1997, A&A, 327, 
966 

Dwek et al 1995, ApJ, 445, 716 

Eaton N., Adams D. J., Giles A. B., 1984, MNRAS, 208, 241 

Epchtein N., 1998, in N. Epchtein ed. The impact of near-infrared 
sky surveys on galactic and extragalactic astronomy, Kluwer 
Academic Publishers, Dordrecht, p. 3 

Freeman K.C., 1987, ARA&A, 25, 603 

Freudenreich H. T., 1998, ApJ, 492, 495 

Garzon F., Hammersley P. L., Mahoney T., Calbet X., Selby 
M. J., Hepburn 1., 1993, MNRAS, 264, 773 

Garzon F., Lopez-Corredoira M., Hammersley P. L., Mahoney T., 
Calbet X., Beckman J. E., 1997, ApJ, 491L, 31. 

Gilmore G., Wyse R. F. G., Kuijken K., 1989, ARA&A, 27 ,555 

Hammersley P. L., Garzon F., Mahoney T., Calbet X., 1994, MN- 
RAS, 269, 753 

Hammersley P. L., Garzon F., Mahoney T., Calbet X., 1995, MN- 
RAS, 273, 206 



© 0000 RAS, MNRAS 000, 000-000 



14 P. L. Hammersley et al. 



Hammersley P. L., Garzon F., Mahoncy T., Calbet X, Lopez- 
Corredoira M., 1996 in D.L. Block and J.M.Grccnbcrg cd., 
New Extragalactic Perspectives in the New South Africa, 
Kluwer Academic Publishers, Dordrecht, p. 534 

Hammersley P. L., Garzon P., Mahoney T., Lopez-Corredoira M., 
1998, in N. Epchtein ed. The impact of near-infrared sky sur- 
veys on galactic and extragalactic astronomy, Kluwer Aca- 
demic Publishers, Dordrecht, p. 63 

Haywood, M., Robin, A. C, Creze, M. 1997, AA, 320, 428 

Kawara K., Kozasa T., Sato S., Kobayashi, Y., Okuda H., Jugaku 
J., 1982, PASJ, 34, 389 

Lopez-Corredoira M., Garzon P., Hammersley P. L., Mahoney T., 
Calbet X., 1997, MNRAS, 292, 15L 

Lopez-Corredoira M., Garzon P., Hammersley P. L., Mahoney T., 
Calbet X., 1999, MNRAS, submitted. 

Mahoney, T., 1998, in preparation 

Majewski S. R., 1993, ARA&A, 31, 575 

Mikami T., Ishida K., Hamajima K., Kawara K., 1982, PASJ, 42 
223 

Neugebauer C, Leighton R. B., 1969, Two Micron Sky Survey. 
NASA SP-3047, GPO, Washington DC 

Paul E. R., 1993, The Milky Way Galaxy and Statistical Cos- 
mology, 1890-1924. Cambridge University Press, Cambridge. 

Reid N., 1993, in S. R. Majewski, ed. Galaxy Evolution: The 
Milky Way Perspective. San Francisco, ASP, p. 37 

Robin A. C, Creze M., Mohan V. 1992, A&A, 265, 32 

Ruphy, S. et al. 1997, A&A, 326, 597 

Schmidt-Kaler T., 1977, Vistas Astron., 19, 69 

Smith, L. F. & Maeder, A. 1991, A&A, 241, 77 

Skrutskie M., 1998, in N. Epchtein ed. The impact of near-infrared 
sky surveys on galactic and extragalactic astronomy, Kluwer 
Academic Phblishers, Dordrecht, p. 11 

Wainscoat R. J., Cohen M., Volk K., Walker H. J., Schwartz D. E., 
1992, ApJS, 83, 111 [WCVWS] 

Wainscoat R. J., Freeman K. C, Hyland A. R., 1989, ApJ, 337, 
163 

Walker H. J., Cohen M.,Volk, K., Wainscoat R. J., Schwartz D. E., 
1989, AJ, 98, 2163 

Weaver H. 1974, in Highlights of Astronomy, cd G. Contopoulos, 
IAU,p423 

van der Kruit P. C, Scarlc L., 1982, A&A, 110, 79 

van de Kruit P. C, 1988, A&A, 192, 117 

Zhao, H. S, 1997 astrphy9705046 



Figures Bl to B9 show in the left column the cumulative 
star counts from the TMGS for each magnitude between 
mK~5 and 9 plotted as the solid line. The longitude at which 
the strip crosses the plane is given in the upper box. The 
long dashed line shows the total prediction from SKY and 
the dotted line the counts from the SKY disc alone. In the 
right column is shown the ratio of TMGS counts to SKY 
total counts for each magnitude. 

Figures CI to C6 show the cumulative counts after sub- 
tracting the SKY disc from each magnitude between mx =5 
and 9. The longitude at which the strip crosses the plane 
is given in the upper box. The solid line is TMGS - SKY 
disc. The long dashed line is SKY total - SKY disc and the 
dotted line is the SKY bulge. 

Figures Dl to D9 (33 plots) illustrate our interpreta- 
tions of the counts by direct comparisons with the predic- 
tions of SKY every half a degree in latiude along the full 
TMGS scans. Differential counts were scaled per square de- 
gree per magnitude, using entirely independent TMGS data 
in 0.2-mag wide bins. The counts arising from the separate 
geometric components are distinguished as follows: solid line 

- total count; dotted or less heavy line - disc; long dashes 

- spiral arms, local spurs and Gould's Belt; short dashes - 
bulge; short dashes and dots - molecular ring; long dashes 
and dots - halo. The position for each area is given in each 
plot. See section 3 and Tables 2 and 3 for more details. Due 
to the number of plots these are only available electronically. 



APPENDIX A: THE PLOTS SHOWING THE 
COMPARISON BETWEEN TMGS AND SKY 

The appendix presents the majority of the the plots showing 
the detailed comparison between TMGS and SKY. From 
these plots it is possible to see how the relative importance 
of the various Galactic components vary with magnitude 
and position. 

Figures Al to A7 (9 plots) illustrate our interpretations 
of the counts by direct comparisons with the predictions of 
SKY. The plots are for a positions every few degrees along 
each scan (the full data is presented in figure D). Differential 
counts were scaled per square degree per magnitude, using 
entirely independent TMGS data in 0.2-mag wide bins. The 
counts arising from the separate geometric components are 
distinguished as follows: solid line - total count; dotted or 
less heavy line - disc; long dashes - spiral arms, local spurs 
and Gould's Belt; short dashes - bulge; short dashes and 
dots - molecular ring; long dashes and dots - halo. The 
position for each area is given in each plot. See section 3 
and Tables 2 and 3 for more details. 
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Figure Ala: TMGS differential star counts (diamonds) plotted on top of the SKY predictions 

for the strip crossing the plane at I = —1°. 
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Figure Alb: TMGS differential star counts (diamonds) plotted on top of the SKY predictions 

for the strip crossing the plane at I = —1°. 
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Figure A2a: TMGS differential star counts (diamonds) plotted on top of the SKY predictions 

for the strip crossing the plane at / = 7°. 
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Figure A2b: TMGS differential star counts (diamonds) plotted on top of the SKY predictions 

for the strip crossing the plane at / = 7°. 






1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 ♦ 1 y 


1 


: 




yf 


: 


















; 


CO 


/' * 


'-. 


- 




*/ 


/- 


: 


J=! ^ 


/ 


: 






♦ /* 




: 




♦/* '' 








/. 








/• 










/ - 


; 
























: 




i J / 


: 












// 






: 


/ 


'-'* /// 


: 




// 






















- 




, , , i/x , , , 1 , , , 


- 









1 1 1 1 ,| 1 1 1 1 l/i /^ 






o 
o 


«3 
O 




-_ 




o 


O 


* / / ■■' 


- 








' //''' 


^- 














II 




♦ / .'■' 
















J2 


-- 


• / /' / x'' 


- 




























.* /' -'■ 


- 


L 






*// / X 


Jl 




















/ ^ 










/ /- 




^ 


// 


^-^ 


/ - 




♦ 










♦, 














f '/ 












_ 










- 




// 




- 








/ ''''■*■ 




















/ ^''' 


- 




1 .'i 


1 1 


\ J\ \ \ \ \ \/a \ \ \ \ 


1 1 









, 1 11 11 1 11 ,», 1 J 






o 




, 






g; ♦ 








'^ ♦ / 








o // * ; 


- 






cv? y* 








♦//*' /' 


. 












" 




\ 




^ 


— * /''' 


J - 










- 




♦ //' 


y - 






































*/''' 










// - 














♦ /.' / / 








W ^ 


/" 




















// 


♦ / y // 


: 


V 


// 




J 














/ ^ •'' 


" 






, 1 , / , , 1 ,//', , 1 , 


, , , " 





' 1 


' ' 


X - 


o 


Ln 




/ : 


o 


n 




♦ /' 


ID 


f^ 




//' ♦ ,' - 








/ ./j 








*// 


^ 


_ 




// * 








V/ 








//■ 1 / / 










- 






X '''' ■' '-'''''" 










- 




J 


♦ / ;'' / /'' : 






5" 








♦ 










/'■ /" 
















* ,v^ 








/ ♦ /// 
























■ '■(/ - 




, -r' 




y:::/ ,,,,^ 






Qfl 






^ 



(2_§9p ^_Si3m) N Sot 



Figure A3: TMGS differential star counts (diamonds) plotted on top of the SKY predictions for 

the strip crossing the plane at I = 21°. 
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Figure A4: TMGS differential star counts (diamonds) plotted on top of the SKY predictions for 

the strip crossing the plane at I = 27°. 
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Figure A5: TMGS differential star counts (diamonds) plotted on top of the SKY predictions for 

the strip crossing the plane at I = 31°. 
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Figure A6: TMGS differential star counts (diamonds) plotted on top of the SKY predictions for 
the strips crossing the plane at I = 37° and I = 53°. 
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Figure A7: TMGS differential star counts (diamonds) plotted on top of the SKY predictions for 

the strip crossing the plane at I = 65°. 
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Figure Bl: Right: cumulative star counts; TMGS (solid line); SKY total (dashed); SKY disc 

(dotted). Left: Ratio of TMGS to SKY total 
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Figure B2: Right: cumulative star counts; TMGS (solid line); SKY total (dashed); SKY disc 

(dotted). Left: Ratio of TMGS to SKY total 
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Figure B3: Right: cumulative star counts; TMGS (solid line); SKY total (dashed); SKY disc 

(dotted). Left: Ratio of TMGS to SKY total 
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Figure B4: Right: cumulative star counts; TMGS (solid line); SKY total (dashed); SKY disc 

(dotted). Left: Ratio of TMGS to SKY total 
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Figure B5: Right: cumulative star counts; TMGS (solid line); SKY total (dashed); SKY disc 

(dotted). Left: Ratio of TMGS to SKY total 



1000 



500 - 



500 
400 
300 
200 
100 

200 
»150 

§100 

o 

" 50 

50 
40 
30 
20 
10 

12 
10 
8 
6 
4 
2 




1 — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



1 

mk = 9 




T 

/=37' 



I I I I I I . I I I I I I I I : 

mk = 8 




n \ \ I \ \ r 



10 10 

latitude 



2.5 ^ 



1.5 r 

1 r 

0.5 \ 

2.5 ~\ 

2 \ 

1.5 ~l 

1 \ 
0.5 ~\ 

2.5 r 

O ^ : 

0.5 I 

2.5 r 

2 r 

1.5 ^ 

1 \ 

0.5 ^ 

2.5 r 

2 ^ 

1.5 r 
1 \ 

0.5 ^ 



mk = 9 



1=37' 




mk = 8 



I I I 
mk = 7 



mk = 6 



mk = 



I I I I I i : 




I I L 




I I I 1 r M I i : 




I I I I 



10 10 

latitude 



Figure B6: Right: cumulative star counts; TMGS (solid line); SKY total (dashed); SKY disc 

(dotted). Left: Ratio of TMGS to SKY total 
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Figure B7 Right: cumulative star counts; TMGS (solid line); SKY total (dashed); SKY disc 
(dotted). Left: Ratio of TMGS to SKY total 
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Figure B8 Right: cumulative star counts; TMGS (solid line); SKY total (dashed); SKY disc 
(dotted). Left: Ratio of TMGS to SKY total 
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Figure B9 Right: cumulative star counts; TMGS (solid line); SKY total (dashed); SKY disc 
(dotted). Left: Ratio of TMGS to SKY total 
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Figure C1& C2 TMGS - SKY disc (solid); SKY total - SKY disc (dashed); SKY bulge (dotted) 
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Figure C3 & C4 TMGS - SKY disc (solid); SKY total - SKY disc (dashed); SKY bulge (dotted) 
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Figure C5 & C6 TMGS - SKY disc (solid); SKY total - SKY disc (dashed); SKY bulge (dotted) 



